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ABSTRACT 

The temperature field produced by a cryosurgical probe 
embedded in a clear gelatin-water test medium has been 
experimentally investigated. The temperature field emanating 
from the cryoprobe was measured using a thin mylar sheet 
sprayed with liquid crystals, a temperature sensitive materi- 
al. Probe tip temperatures varying from -36 °C to -117 ®C 
were studied. The experimental results compared within 
approximately 9% of a one dimensional analytical solution 
for predicting ice growth rate. 
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NOMENCLATURE 
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specific heat in the unfrozen region cal/g/®C 
specific heat in the frozen region cal/g/°C 
thermal conductivity of the unfrozen region 
cal/sec/cm/®C 

thermal conductivity of the frozen region 
cal/sec/cm/®C 

latent heat of fusion cal/g 
radial distance cm 
radius of cryoprobe cm 

radial location of the frozen/unfrozen inter- 
face cm 
time sec 

temperature in the unfrozen region ®C 
temperature in the frozen region ®C 
initial medium temperature ®C 
phase change temperature ®C 
cryoprobe tip temperature ®C 
density g/cm^ 
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I. 



INTRODUCTION 



The use of low temperatures as a medical tool is not 
new. As early as 3500 B.C. the Egyptians used ice as a local 
anesthetic by applying it to wounds and other ailments. 
However, it was not until 1851 that extreme cold was used as 
a surgical procedure. Dr. James Arnott, a British physician, 
used a small container filled with an ice-brine mixture at 
-24®C as a surgical instrument. The container was placed in 
direct contact with skin cancers and was found helpful in re- 
stricting the growth of the cancerous tissue. In 1961, with 
the development of the first cryosurgical probe for the treat 
ment of Parkinson's Disease, cryosurgery became an establish- 
ed surgical procedure. Dr, Irving Cooper of the Saint Barna- 
bas Hospital in New York is acknowledged as the founder of 
the present day cryosurgical technique. Cooper developed a 
small multi-tubed probe with a hemispherical tip. The tem- 
perature of the tip could be controlled from body temperature 
37®C, to -196®C by circulating liquid nitrogen through it 
HI . 

The cryoprobe has primary and secondary applications in 
treating diseases of the brain. The major use has been in 
the production of discrete destructive thalamic lesions in 
the surgical treatment of Parkinson's disease and related 
disorders. Unfortunately, the application of the cryosur- 
gical technique to the management of benign or malignant 
tumors has been severely limited due to an inability to 
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predict lesion growth rate and ultimate lesion size [17] , 
There have been several analytical models developed for 
predicting lesion growth rate and size [1, 6, 7]. However, 
these models have yet to be verified experimentally. 

The objectives of this thesis were: 1) to experimen- 

tally determine the growth rate of the frozen region around 
a cryoprobe 2) to investigate the feasibility of using 
liquid crystals, a material that indicates temperature 
through a change in color, as a temperature sensing device 
in the unfrozen region 3) to compare the experimentally 
determined temperatures and ice growth rates with existing 
analytical predictions. 

Two cryoprobes, 0.5 inches and 0.187 inches in diameter, 
were designed and built along the lines of the standard 
Linde cryoprobe. The larger probe, which was intended as a 
prototype, produced ice regions that were approximately 
spherical in nature and easily studied. The smaller probe, 
due to faulty design, produced irregular regions that did not 
submit to detailed study. A 0.15% gelatin 99.85% water test 
medium was used to simulate tissue. A range of probe tip 
temperatures varying from -36®C to -117 ®C were used to 
create the frozen regions. In all cases studied, experimen- 
tal and analytical predictions of ice growth rates compared 
within the estimated experimental uncertainty of 9%. The 
details of the magnitude and extent of the temperature field 
in the unfrozen region were successfully obtained with the 
liquid crystal material. 
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II. DESIGN OF THE SURGICAL CRYOPROBE 



The general surgical cryoprobe used for necrotizing 
brain tumors consists of three concentric tubes [17] . The 
coolant, which is normally liquid nitrogen (LN 2 ) , flows 
down the inner tube and absorbs heat from the tip, causing 
the coolant to vaporize. The vapor flows back through the 
middle annulus. The outer annulus is evacuated to minimize 
heat transfer between the gelatin and the returning vapor 
(Figure 1) . Probe tip temperature is controlled by the 
coolant flow rate, which is in turn controlled by regulating 
the outlet pressure from the coolant supply reservoir. 

The original desire was to determine the ice formation 
rate on an existing Linde cryosurgical unit. Attempts to 
secure such a unit were unsuccessful resulting in a need to 
construct a unit. Using the basic design of a Linde cryo- 
probe [17] as a model, two probes were designed and fabricated 
from stainless steel hypodermic tubing. The first cryoprobe 
built was a prototype. It was 0.5 inches in diameter and 
5.75 inches in length (Figure 2, 3) . The prototype was 
tested and found to work satisfactorily. The same design 
was applied to a smaller cryoprobe, 0.187 inches in diameter. 
This cryoprobe is the same size as the Linde Cryoprobe PR- 2 
used for brain tumors. When the smaller cryoprobe was test- 
ed, it failed to form a spherical ice ball. instead, an 
elliptical ice ball formed indicating two problem areas. 

The vacuum volume was too small and the tip was not a 
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hemisphere but was shaped more like a cylinder. A new 
design should employ a smaller inner tube for the LN 2 to 
flow through and the tip should be a hemisphere as in the 
case of the large probe. The remainder of this thesis will 
be concerned with the 0.5 inch cryoprobe. 
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FIGURE 1 
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0.5 X 5.75 INCH CRYOPROBE 



FIGURE 2 
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0.187 X 6.56 INCH CRYOPROBE 



FIGURE 3 
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III. EXPERIMENTAL APPARATUS 



The experimental set-up used in the present studies 
was similar to the set-up used by Groff in his studies of 
the resistively heated surgical probe [14] . Briefly, the 
set-up consisted of a 8 inch cubical test cell, with the 
cryoprobe horizontally inserted through a small hole in 
the center of one side of the test cell. The liquid crystal 
material, which was sprayed on a mylar sheet, was placed in 
the horizontal midplane of the cryoprobe. The probe had 
thermocouples welded on the tip and on the midpoint of the 
stem. The test cell was filled with a clear 0.15% gelatin 
/99.85% water test medium, which had the following thermal 
properties 

a. = .0050 cal/sec/cm/®C 

b. k = .0014 cal/sec/cm/®C 

c. L = 79.71 cal/g 

The cryoprobe was connected to a LN 2 source. In order to 
record the visual data presented by the liquid crystals, 
a Graflex Camera was placed directly above the cryoprobe 
tip (Figure 4). 

A. LIQUID NITROGEN SYSTEM 

The LN 2 was stored in a Union Carbide CE-8, which is a 
dewar that uses the pressure of the vaporized LN 2 as the 
driving potential to force the LN 2 through a regulating 
valve. The unit is also equipped with a heater which can 
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EXPERIMENTAL SET-UP 



be used to vaporize the LN 2 to increase the pressure as 
necessary. The maximum working pressure of the system is 
10.5 psig. The regulating valve is a needle valve that 
allows the operator to regulate the flow of the LN 2 (Figure 
5) . 

To transfer the LN 2 from the dewar to the cryoprobe a 
connecting tube was designed and built (Figure 6) . The tube 
is essentially two concentric tubes. The LN 2 flows through 
the inner annulus and the outer annulus is evacuated. The 
entire tube is encased in neoprene insulation. 

B. LIQUID CRYSTALS 

There are three types of liquid crystals; smectic, 
nematic and cholesteric [12, 13, 15, 16] . This experimental 
study used only the cholesteric type. A cholesteric liquid 
crystal responds to changes in temperature by sequentially 
passing through the complete visual spectrum (red through 
violet) in fractions or multi-degrees, depending on which 
cholesterol esters comprised the formulation. This color 
phenomenon is reversible. A very important point is that 
at a certain temperature a given material or combination of 
materials will always exhibit the same color. The light 
scattered by the liquid crystals represents only a fraction 
of the incident light. The remaining portion of the incident 
light is transmitted by the liquid crystals. Therefore, an 
absorptive black background must be used to prevent reflec- 
tion of the transmitted light, thereby enhancing the color 
resolution. 
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FIGURE 6 
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since the colors exhibited by the cholesteric liquid 
crystal are unique for a specific temperature, the quanti- 
tative measurement of temperature is possible to an accuracy 
of approximately 0.1®C, while the speed of response of the 
liquid crystals is less than 0.2 seconds. 

C. DESCRIPTION OF LIQUID CRYSTAL SHEETS 

Mylar sheeting was chosen for the foundation of the 
liquid crystal, because it is thin, pliable and has a small 
thermal conductivity, which for this system results in 
negligible heat transfer along the sheet. 

The mylar sheet was cut into a 7 15/16" square with a 
4 1/8" X 1/2" notch cut into the sheet for the experimental 
runs (Figure 7) . The liquid crystals were applied to the 
mylar sheets with a Paasche, Model H "3 in 1" air brush. 

The steps taken in preparing the liquid crystal sheets were: 

A. Thoroughly clean both sides of mylar sheet with 
acetone. 

B. Use liquid slurry received from NCR without dilution. 

C. Spray first coat of liquid crystals on one side of 
the clean mylar surface allowing the crystal to 
dry for 40 to 60 minutes. 

D. Spray second coat and allow to dry for 20 minutes. 

E. Spray third coat and allow to dry for one hour. 

F. Apply the first coat of Testor*s flat black paint 
on the dryed liquid crystals and allow to dry for 
30 minutes. 
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LIQUIL' CRYSTALS SPRAYED ON A MYLAR 
SHEET, THE SHEET IS ENCASED IN THE 
I.IQUID CRYSTAL FRAME. 



FIGURE 7 
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G. Apply the second coat and let dry for one hour. 

H. Apply two coats of Rex Polyurethane Interior- 
Exterior 77-5 over the black paint paying particular 
attention to the edges. This will water proof the 
liquid crystal sheet. 

The mylar sheet for the experimental run was divided 
into five areas (Figure 8) . The W-5, W-10, R-15 and R-19 
areas were covered with masking tape. W-0 crystals were 
applied as described earlier, but only one coat of Testor's 
flat black paint was used. The W-5 area was uncovered next 
and the V7-5 crystals were sprayed. Then one coat Of Testor’s 
flat black paint was applied to the W-5 area. The process 
was repeated for the W-10, R-15 and R-19 areas. A second 
coat of Testor's flat black paint v/as then applied to the 
entire sheet and allowed to dry. The sheet was then water 
proofed. The reader is referred to Appendix A for calibra- 
tion technique and data. 

D. TEST CELL 

The test cell was composed of two parts. The first was 
the liquid crystal sheet frame. The frame, which supported 
the liquid crystal sheet, was 7 15/16 inches square and was 
made up of two plexiglass squares that were screwed together 
with the liquid crystal sheet in the middle (Figure 8) . On 
one side of the frame was a 0.5 inch diameter hole through 
which the cryoprobe fit. 

The second part was a clear plexiglass 8 inch cube. 

Four inches down from the top of the cube was a 0.5 inch 
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wide shelf that supported the liquid crystal frame (Figure 9) , 
A 0.5 inch diameter hole in the side of the box corresponded 
to the hole in the frame. When the cryoprobe was placed 
through the hole in the box and the frame, the longitudinal 
axis of the cryoprobe was in the same plane as the liquid 
crystal sheet. 
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EIGHT INCH CUBICAL TEST CELL 



FIGURE 9 
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IV. EXPERIMENTAL PROCEDURE 



The liquid crystal frame was placed in the box. The 
cryoprobe was pushed through the hole in the box and the 
frame. A copper/constantan thermocouple was welded on the 
tip of the cryoprobe. To record the temperatures, a Honey- 
well Electronik 194 dual strip recorder with ice junction 
was used. 

The gelatin was made and allowed to cool for six hours 
(Appendix C) . After six hours the gelatin was cool enough 
to pour into the assembled box without thermally shocking 
the liquid crystals. The gelatin then was allowed to come 
to room temperature. The LN 2 dewar was connected to the 
cryoprobe by means of an insulated double wall tube (Figure 
6) . A vacuum pump was connected to the cryoprobe and the 
double wall tubing. A pressure of approximately 10mm of 
mercury was maintained in both probe and tubing. The Graf lex 
camera (Appendix B) was positioned directly over the cryo- 
probe and all necessary adjustments were made. The ther- 
mocouples were connected to the strip recorder. The pressure 
on the LN 2 dewar was between 3 and 9.5 psi. 

At time equal zero the LN 2 was turned on. The tip 
temperature was monitored with the recorder. When the tip 
temperature reached a predetermined temperature, the controls 
on the dewar were adjusted so the tip temperature remained 
constant. Pictures were taken at times of 0, 1, 3, 5,. 10, 

15, 30 and every 15 minutes thereafter until the 21.5®C 
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isotherm moved to the edge of the box. At this point in 
time the experiment was terminated. 

The photographs recorded the isotherms and ice ball 
diameter. Using photo interpolation (knowing the probe 
diameter to be 0.50 inches) the actual ice ball diameter 
and the distance from the cryoprobe tip center to the 
various isotherms were determined. 

Six experimental runs were made varying the probe tip 
temperatures as shown below: 

Experiment Tip Temperature Initial Gelatin 

Temperature 



1. 


-36®C 


22®C 


2. 


-53®C 


21®C 


3. 


-65°C 


21°C 


4. 


-15^C 


20°C 


5. 


-97®C 


25.5®C 


6. 


-117*^0 


22.5°C 



In all cases the phase change temperature was -2.0®C. 
This was found by allowing the strip chart recorder to run 
for one hour after the experiment ended. The probe tip 
thermocouple steadied on -2®C and remained there until the 
ice ball had completely melted away from the cryoprobe. 

Figure 10 shows a segment of chart paper for a typical 
experimental run. The first half is the response when the 
LN 2 is turned on and the second half is the response when 
the LN 2 is turned off. 

Figures 11, 12, 13, and 14 show the ice ball forming in 

the gelatin with the liquid crystals visually presenting 
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the isotherms. The relationship between the colors and 



temperatures are: 
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0.5 INCH CRYOPROBE EMBEDDED IN THE 
GELATIN SURROUNDED BY THE LIQUID CRYSTAL 
SHEET AT TIME EQUAL ZERO. SEE FIGURE 7 
FOR COLOR TEMPERATURE RELATIONSHIPS. 



FIGURE 11 
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0.5 INCH CRYOPROBE EMBEDDED IN THE 
GELATIN SURROUNDED BY THE LIQUID 
CRYSTAL SHEET AT TIME EQUAL FIVE MINUTES. 

SEE FIGURE 7 FOR COLOR TEMPERATURE RELATIONSHIPS. 



FIGURE 12 
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0.5 INCH CRYOPROBE EMBEDDED IN THE GELATIN ♦ 

SURROUNDED BY THE LIQUID CRYSTAL 
SHEET AT TIME EQUAL THIRTY MINUTES 
SEE FIGURE 7 FOR COLOR TEMPERATURE RELATIONSHIPS. 



FIGURE 13 



33 




0.5 INCH CRYOPROBE EMBEDDED IN THE GELATIN 
SURROUNDED BY THE LIQUID CRYSTAL SHEET 
AT TIME EQUAL SIXTY MINUTES. SEE FIGURE 7 
FOR COLOR TEMPERATURE RELATIONSHIPS. 



FIGURE 14 
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V. THEORETICAL ANALYSIS 



When a cryoprobe is inserted into tissue and brought to 
a cryogenic temperature, a thermal field develops in the 
tissue. A lesion is formed as the tissue in the immediate 
vicinity of the cryoprobe undergoes a phase change. The 
freezing front will propagate through the tissue until steady 
state is reached within the tissue. The specific heat, den- 
sity, thermal conductivity and latent heat of fusion, as well 
as blood flow and metabolism all contribute in controlling 
ultimate lesion size and growth rate for a given cryoprobe 
diameter and particular tip temperature. Therefore, to 
predict lesion size in a structure such as a human brain, 
it is necessary to know not only the cryoprobe tip temperature, 
but also the thermal properties of the structure immediately 
adjacent to the cryoprobe tip [4] . 

There have been several models [1, 6, 7] that describe 
the transient behavior of a cryolesion. Cooper and Trezek 
[6, 7] present an approximate solution to the problem of 
transient freezing around cylinders and spheres embedded in 
tissue. Their model included terms for the thermal proper- 
ties, as well as for blood flow and metabolism. While the 
blood flow and metabolism terms are important, they have 
been disregarded in this study. In order to take the 
phenomenon into account, a live brain is necessary. Since 
the school lacked the proper facilities for operating on 
experimental animals, it was decided to experimentally 
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investigate the rate of ice formation produced by a cryo- 
probe embedded in gelatin. The experimental results were 
compared to the analytical solution of Cooper and Trezek 
to ascertain if their basic premise (neglecting heat capacity 
effects in both the frozen and unfrozen phases) is valid. 

The following is a brief development of the solution 
presented by Cooper and Trezek without metabolism and blood 
flow terms. 

The energy equation (see Nomenclature for definition of 
terms) 

pc II = kV^T (1) 



is used to describe the temperature field around the cryo- 
probe. The energy equation will take the following form 
when used to describe one dimensional heat transfer in the 
frozen and unfrozen regions surrounding a spherical probe. 

Unfrozen region 



pc 



at ^2 3ry 

A. (r2 

3r V 3r/ 



Frozen region 

3T _ „2 






with the following boundary conditions 



(2) 



( 3 ) 



r 



= 



( 4 ) 



r = 

r oo 




( 5 ) 

( 6 ) 



r = 





+ pL 



dt 



( 7 ) 



36 



By introducing the following non-dimensional groups, 



T- - T 
= .J o 

f - T 
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the normalized energy equation becomes: 
frozen region 
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r2 3R 
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unfrozen region 
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(9) 



with the boundary conditions 4, 5, 6, and 7 becoming 



R = 


1 
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(10) 
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r* 
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R = 
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36 , 
3 ^ 
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(13) 



Analytical solutions of equation (8) and (9) satisfying 
conditions (10) , (11) , (12) , and (13) are not known. However, 

if we assume 



k c . (T - T ) c (T - T ) 

and — = are small, that 



k^L 

is, if the heat capacity effects are small compared to latent 
heat effects (8) and (9) take the approximate form: 
frozen region 

90, 



1 

r2 



R2 

9R \ 3R 



= 0 



(14) 
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unfrozen region 



8 

9R 






(15) 



the solutions to equations (14) and (15) are 
unfrozen region 

e * 

6 = . 



frozen region 






1 + 






pc 






(16) 



(17) 



Equations (18) and (19) are substituted into equation (13) 
to find V 



dr* 

dx 






1 -e 



pc 



\r + (r*-ll 



Further rearrangement of equation (18) yields 



dx = 



_ r* (1 - r*) 



(r*+n ) 



dr* 



(18) 



(19) 



where 



f 

"F 



' pc ’ 



dr * 

at steady state - 0, this necessitates from equation 

(19) that 

r* = - n (20) 

at steady state. Integrating (19) between the limits 1 and 
-n yields: 



X = (n+1) (r*-l) + n(Ti+l)ln 



n + 1 



(r*-l) (21) 



Ti+r* 2 

Using equation (20) and (21) an analytical prediction 
of the rate of growth of the ice ball can be made. A 
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FORTRAN program (Appendix E) was written that uses equation 
(20) and (21) to predict the rate of lesion growth. The 
program also predicts time to steady state and maximum 
lesion size for the particular cryoprobe size and tip tem- 
perature . 

When using the FORTRAN program written in conjunction 
with ..this thesis the actual radius of the hemispherical 
cryoprobe tip is used vice a modified radius. At the 
beginning of this study it was felt that perhaps the radius 
should be the radius of an imaginary sphere that has the 
same surface area as the hemispherical cryoprobe tip since 
the theoretical analysis assumed a spherical case. However, 
the best agreement between the theoretical solution and 
experimental data occurred when the actual hemispherical 
radius of the cryoprobe tip was used in the FORTRAN program. 
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VI . RESULTS 



Correlation between the experimental data and the theory 
was excellent. The theoretical curves for the six different 
cryoprobe tip temperatures, as calculated by the FORTRAN 
program (Appendix E) , are illustrated in Figures 15, 16, 17, 
18, 19, and 20. The experimental data is superimposed on 
the graphs to demonstrate the agreement. For the first 
hour the correlation is almost exact, then the experimental 
data starts to rise above the theoretical curve, but remains 
within the uncertainty bounds (Appendix D) . The rise in the 
experimental data can be explained by the theoretical 
boundary conditions not being met after the first hour. 

That is, the 19®C isotherm had moved to the edge of the test 
cell and the gelatin no longer qualified as an "infinite" 
medium. 

Figure 21 is a typical graph of the nondimen sional 
time vs. nondimensional radius for a tip temperature of 
-25®C, The solid line represents the theoretical solution 
and the dashed lines represent the uncertainty bounds. 

The reader is referred to Appendix D for details on the uncer- 
tainty analysis. 

The transient temperature field shows the details of 
a lesion formation. An example is illustrated in Figure 
22. The discontinuity at -2®C, the phase change temperature, 
indicates the lesion size at a particular time after initiat- 
ing cooling. The experimental points on the graph were 
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COMPARISON OF THEORETICAL AND EXPERIMENTAL 
PREDICTION OF ICE FRONT LOCATION, r* , AS A 
FUNCTION OF TIME, t, FOR A PROBE TEMPERATURE 
OF -36°C. 

FIGURE 15 
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COMPARISON OF THEORETICAL AND EXPERIMENTAL PRE- 
DICTION OF ICE FRONT LOCATION, r*, AS A FUNCTION 
OF TIME, T, FOR A PROBE TEMPERATURE OF -53 ®C. 



FIGURE 16 
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COMPARISON OF THEORETICAL AND EXPERIMENTAL PRE- 
DICTION OF ICE FRONT LOCATION, r*, AS A FUNCTION 
OF TIME,!, FOR A PROBE TEMPERATURE OF -65®C. 



FIGURE 17 
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COMPARISON OF THEORETICAL AND EXPERIMENTAL PRE- 
DICTION OF ICE FRONT LOCATION, r*, AS A FUNCTION 
OF TIME, T, FOR A PROBE TEMPERATURE OF -75°C. 



FIGURE 18 
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COMPARISON OF THEORETICAL AND EXPERIMENTAL PRE- 
DICTION OF ICE FRONT LOCATION, r*, AS A FUNCTION 
OF TIME, T, FOR A PROBE TEMPERATURE OF -97®C. 



FIGURE 19 
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COMPARISON OF THEORETICAL AND EXPERIMENTAL PRE- 
DICTION OF ICE FRONT LOCATION, r*, AS A FUNCTION 
OF TIME, T, FOR A PROBE TEMPERATURE OF -177 °C. 



FIGURE 20 
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GRAPH OF THE UNCERTAINTY BOUNDS OF THE THEORET- 
ICAL CURVE AND UNCERTAINTY LIMITS OF THE EXPERI- 
MENTALLY DETERMINED DATA AT A CRYOPROBE TIP 
TEMPERATURE. OF -75®C. 

FIGURE 21 
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GRAPH OF THE TRANSIENT TEMPERATURE FIELD AT 
TIMES OF 5, 30, 90, 210 MINUTES AT A CRYOPROBE 
TIP TEMPERATURE OF -75°C. 



FIGURE 22 
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determined from the liquid crystal sheet. The experimentally 
determined points show that the temperature gradient in the 
unfrozen phase is much sharper than the theoretically pre- 
dicted temperature gradient. This difference is attributed 
to the heat capacity effect which is neglected in the theo- 
retical analysis. However, the thermal capacitance does 
not seem to be important in determining the ice growth rate. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



The use of liquid crystals afforded an unique opportu- 
nity to visually witness the rate of movement of various 
isotherms. The analytical prediction of the lesion growth 
rate, as presented by Cooper and Trezek, shows good agree- 
ment with experimentally determined data. The assumption 
made in the theoretical analysis, that the heat capacity 
in the frozen phase is negligible, seems reasonable. How- 
ever, the thermal capacity in the unfrozen region does have 
a definite influence on the extent of the temperature field. 

The ability to control the tip temperature was difficult 
at higher temperatures (-36°C, -53°C) and relatively easy 
at lower temperatures (-117 °C) . It is recommended that an 
automatic control system be developed to control the cryo- 
probe tip temperature. 

The smaller cryoprobe should function properly if the 
tip were redesigned so as to be hemispherical. Also when the 
small cryoprobe is used in an experimental run, low LN 2 
pressure (1.5 to 3.5 psig) should be used. 
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APPENDIX A 



Liquid Crystals 

For this project seven different temperature ranges of 
the cholesteric type were used. The manufacturer was 
National Cash Register (NCR) . The crystals used were: 

(NCR Notation) 



A. 


W-0 


B. 


W-5 


C. 


W-10 


D. 


R-15 


E. 


R-19 


F. 


R-23 


G • 


R-27 



W - wide temperature range -5®C 
R - regular temperature range 4®C 

No information was received as to the exact temperature at 
which color changes would occur nor the exact meaning of 
the ranges 5®/4®C. Thus it was necessary to calibrate the 
crystals prior to any experimentation to find their ranges 
of utilization. 

Two calibration mylar sheets were prepared. Both sheets 
were four inch square. The same technique used for preparing 
the liquid crystal sheet used in the experiment was used for 
preparing the calibration sheets. Each sheet was divided 
into fourths. The first sheet contained the W-0, W-5, 
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W-10, and R-15 crystals. The second sheet was similarly 
divided with R-19, R-23, and R-27 applied. Each sheet was 
calibrated separately in the Rosemount, model 913A, calibra- 
tion bath (constant temperature bath) . The bath was lowered 
to 0°C. The first sheet, containing W-0, W-5, W-10 and 
R-15, was placed in the bath. The bath temperature was 
monitored by a Rosemount Commutating Bridge, model 92A, and 
galvanometer arrangement with a Rosemount Platinum Resistance 
Temperature Standard 162C. Throughout the calibration pro- 
cedure one observer read the color changes while the other 
read the commutating Bridge. The following results were 
obtained . 

Calibration of Liquid Crystals 





Red 


Green 


Blue 


W-0 


1.7°C 


6.2®C 


8.0°C 


w-5 


3.3®C 


9.3°C 


11.0°C 


W-10 


10.0°C 


12.7®C 


15.4®C 


R-15 


15.9“C 


16.5°C 


17.1°C 


R-19 


19.8°C 


U 

o 

00 

o 

CM 


21.5°C 


R-23 


22.4®C 


23.5°C 


24.6®C 


R-27 


26.7«C 


27.9®C 


28.9°C 
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APPENDIX B 



Photographic Techniques 

In order to make a permanent record of the experiments, 
pictures were required. For data taking Polaroid Polacolor 
4x5 land film Packets, color/type 58, were used in a Graf lex 
camera (4x5 press Graf lex) . The experimental runs were 
made at night. The lighting was supplied by a Tungsten 
3200°K light. The film used was a daylight type necessitat- 
ing the use of a Tiffen Photar Filter SOB series 6 for color. 
When taking pictures for publication quality, 35mm high 
speed Ektachrome color film was used in a Honeywell Pentax 
H~3. This film was designed for use with the Tungsten 
3200°K lights. 

In both cases the technique used for taking the pictures 
was the same. The camera, on a tripod, was placed directly 
over the cryoprobe. The camera was focused and settings 
were made. The Tungsten light was placed approximately six 
feet away from the cryoprobe. When a picture was desired, 
the Tungsten 3200®K light was turned on, the room lights 
turned off and the picture taken. As soon as the picture 
was taken, the Tungsten light was turned off. The total 
time the Tungsten light was left on was approximately ten 
seconds. It was felt that the heat generated by the light 
had little or no effect on the freezing rate. 

The general photographic techniques were used when 
taking pictures through the gelatin. The highest f stop 
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with the slowest shutter speed was used to get the greatest 
depth of field. The main drawback with the Polaroid film 
was that it was a daylight film with a reasonable ASA, but 
when the filter was added the ASA was lowered to 24. Color 
reproduction was considered adequate for data taking. 
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APPENDIX C 



Preparing the Gelatin 

To make the gelatin the following recipe was used: 

a. Heat 6800 ml of v/ater to boiling 

b. Soak 120 gms of Schilling *s Gelatin in 1200 

ml of cold water for 3 minutes. Do not disturb 
during this time. 

c. Pour boiling water into the gelatin mixture. 

d. Stir gently while pouring in the hot water. 

If stirred too vigorously a foam will appear on 
the surface. 

e. After all hot water is in gelatin mixture, stir 
to insure all gelatin is thoroughly mixed. 

f. Leave gelatin in mixing container and put 

in refrigerator (initial temperature 32°F) for 
3 hours. 

g. Take out of refrigerator and let stand at room 
temperature for 3 hours. 

h. Pour into experimental box. 

i. Let stand at room temperature for 12 hours. 
Gelatin is now ready for use. 
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APPENDIX D 



The uncertainty in a function f(x,y,z) computed from 
measured values of several independent quantities (x,y,z) 
may be estimated as 




where is the estimated uncertainty in f due to estimated 
uncertainties of in x, y, z respectively [ 2 , 11] . 



Such an uncertainty analysis was carried out on the 
- T^\ k(T^ - T t 

I ^ 



quantities n = 



V 



T - T 
pc o , 



-1 and T = 



_£C_ 



pLr^ 



to ascertain the influence of experimental uncertainties on 
the theoretical and experimental determinations of r* . 
Applying (1) to n and t it is found; 

» ■'/ft / • (» I ■ (5^)' ■ 



and 




The following are the estimated true values and uncer- 
tainties for the experimental run with the probe temperature 



held at -75®C. 

T - T = -77®C ± 1°C 
pc s 



20:1 



T - T = -22®C ± 1°C 20:1 

pc o 
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k = 0.0014 + .0001 

C cm sec 



20:1 



0.0050 ± .0001 



Cal 

C cm sec 



20:1 



p = 1.00 ± .01 22 - 

Cm3 



20:1 



L = 79.71 ± .01 



Cal 

gm 



= 0.025 ± .002 Cm 



20:1 

20:1 



T ~ T 
pc s 

r - T 
pc s 



0.01 



= 0.01 

p 



T - T 
pc o 

‘’t - T 

pc o 



0.05 



= 0.01 



Zk 

k 



0.07 



= 0.01 



"1^ 



= 0.02 



The results of the uncertainty analysis indicated that 
for this particular run: 

0 ) 0 ) 

— =0.09 and -2. =0.09 

T n 

Similar calculations were made on the other five experi- 

U) (1) 

mental runs and in all cases — and — fell between 0.08 and 

T n 

0.09. A typical uncertainty plot is shown in Figure 21. 
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